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ABSTRACT 

Topic of this paper is the time-evolution of intrinsic correlations of galaxy ellipticities due 
to peculiar motion. In our model, the galaxy ellipticities are determined from the angular 
momentum of their host haloes, which can be computed from the fluctuations statistics of a 
Gaussian random field. Subsequent peculiar motion distorts the ellipticity field and causes 
changes in the ellipticity correlations. Using analogies between this problem of shifted el- 
lipticity tensors and the displacements of polarisation tensors in gravitational lensing of the 
cosmic microwave background we compute £^-mode and 5-mode spectra of the time-evolved 
ellipticity field, where the displacements are modelled with first and second order Lagrangian 
perturbation theory. For EUCLID, ellipticity correlations are decreased on large multipoles 
i ^ 1000, amounting to up to 10% in the £^-mode spectrum C|(^) and up to 60% in the B- 
mode spectrum C|(^) at ^ ^ 3000 due to the dispersing eff'ect of peculiar motion. £/5-mode 
conversion in analogy to CMB -lensing is present but small. We conclude that distortions of 
the ellipticity field due to peculiar motion is not aff'ecting the prediction of ellipticity models 
on the scales relevant for lensing in the case of EUCLID 's galaxy distribution, but should 
aff'ect larger scales for surveys at lower redshifts. 
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1 INTRODUCTION 

Weak gravitational lensing by the cosmic-large scale structure is a 
tool for investigating the fluctuations statistics of the cosmic density 
field and its dependence on the und erlying cosmological model (for 
reviews, see Bartelmann & Schne ider 2001 ; Bartelmann 2010). Fu- 
ture weak lensing surveys such as EUCLID, DES, LSST and IDEM 
are designed to yield sub-percent accuracy on the set of cosmo- 
logical parameters in a dark energy cosmology by measuring the 
correlation function of the gravitationally sheared galaxy ellipticity 
field. A common assumption in weak cosmic lensing is that those 
ellipticities are intrinsically uncorrected and that the only corre- 
lating eff'ect is weak lensing, because the light from neighbouring 
galaxies has to transverse the same cosmic tidal fields, leading to a 
correlation in change in shape. 

This assumption, however, is challenged on small scales by 
intrinsic alignment eff'ects (for a re view on ang ular momentum 
models and intrinsic alignments, see ISchafei] [2009). Tidal shear- 
ing models for angular momentum built-up in galactic haloes pre- 
dict correlated angular momenta of neighbouring galaxies. If the 
symmetry axis of the galactic disk is aligned with the angular 
momentum direction of the host halo, neighbouring galaxies are 
viewed under correlated angles of inclination such that their ellip- 
ticities appear correlated. This intrinsic alignment eff'ect is impor- 
tant on small scales because the angular momentum correlation is 
comparatively short ranged: it is predicted to be present on scales 
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of about 1 Mpc//z ([Crittenden et all 1200 ll : iNataraian et"aDl200ll : 
ISchafer&Merkelll201lh . 

Intrinsic alignments of galaxies based on angular mo- 
mentum models are a relatively new topic which will un- 
doubtedly attract much interest in the future as the weak 
lensing data sets provided by large-scale lensing observations 
will at the same time help to scrutinise intrinsic alignment 
models. The theory of angular momentum-induced alignments 
teroft & Metz"ie3l200Ql : ICrittenden etaDl200lL l2002l : iMackev et all 
2002) has been applied to describe con tamination of weak 
lensin g data in the convergence spectrum ( Heavens et al.l |200QI : 
„Hevmans & HeavensI l2003l: iHevmans etaP |2004 l2006h and bis- 
pectrum dSemboloni et al.l l2008h . Difi'erent schemes for remov- 
ing the contamination of intrinsic alignment s have been pro- 
posed , from discarding close galaxy pairs (iKing & Schneidej 
I2003L l2002h to specifically desig ned weighting sch emes for 
nulling out their contribution dJoachimi & Schneidej l2008i) or 
amplifying them relative to the weak lensing induced elliptic- 
ity c orrelations (King & Schneider 2003; "joachimi & Schneidei] 



20101) . Resulting biases on cosmological parameter estimation 
if intrinsic alignment s remain uncorrected have been quantified 
teridle & King' '20071; Ijoachimi & Bridl^ |201QI : iKirk et all |201QI : 
Schneider & Bridle 20l3). 

The wealth of structure in the angular momentum and ellip- 
ticity field and their alignment with large-scale tidal fields ha s at- 
tracted much interest fro m the numerical perspective ( Hahn et al.l 
I2OOI ICodisirZll2012h and suggests the question if large-scale 
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tidal field s can be reconstructed using the ellipticity field as a tracer 
(iLee & Pe n 2000, 2001, 2007). 

On a more fundamental level, the investigation of tidal 
shearing mechanisms in diff'erent orders of perturbation theory 
along with the deformation of forming haloes due to tidal forces 
dCatelan & Theuns 'l996b'; 'Catelan' h99^; 'Catelan et al. 2001; 
Catelan & Theuns 19 96a; Catelan et al. 1995; Catelan & Porciani 
200ll : lLee et al.ll2007h . N^umerical verification of the tidal torquing 



theories of angular momentum build-up has been the topic of a 
number of papers ( Catelan & Porciani 200l1; I Bullock et al ] |200ll : 
IPorciani et al.ll2002al lbl: IHahn et all 12007) who agree that the an- 
gular momentum direction can be well described by tidal torquing 
whereas the amount of angular momentum might be overestimated. 
Because ellipticity alignments are only sensitive to the angular mo- 
mentum direction, these studies provide support for using tidal 
torquing models with this particular application in mind. 

The way in which the orientation of a galactic disk is linked to 
the angular momentum direction of the host halo is not very clearly 
cut. In a small number of observations iBailin et all (l2005h found 
mismatches and suggest that direct linking of the symmetry axis of 
the disk to the host halo angular momentum would lead to overesti- 
mation of the ellipticity alignments. This eff'ect is partially covered 
by parameterisation, but unless the relation is better understood, 
angular momentum-based alignment models proved upper limits 
on ellipticity correlations. By now, intrinsic alignments have been 
measured in a number of data sets and have been found at the ex- 
pected levels ("Pen et al. ll200Ql : lLee & Penll2002l : lMandelbaum et all 
I2OO6 ; Hi rata et al. 2007 ), alth ough some studies doubt these claims 
(e.g. lAndrae & Jahnkd l2oTT). 

The point which motivated this paper is the comparatively 
short-ranged correlation of angular momenta and consequently of 
the galaxy ellipticities, which reaches out to distances of about 
1 Mpc/h. If the ellipticity field is this short-ranged, and if it is 
distorted by the peculiar motion of galaxies, by how much do the 
correlations change and on what scales? We will investigate this 
question by employing a formalism based on lensing of the CMB- 
polarisation, by describing peculiar motion with Lagrangian per- 
turbation theory, and consider the EUCLID galaxy sample as an 
application: After a summary of cosmology, structure formation, 
Lagrangian perturbation theory, angular momentum models and el- 
lipticity correlations in Sect. [21 we describe our formalism and the 
results in Sect. [3] Our main findings are summarised in Sect. ID 

As reference model we chose a spatially flat wCDM model 
with Gaussian adiabatic initial perturbations in the cold dark matter 
distribution. Specifically, parameters were chosen to be = 0.25, 
n, = 1,0-8 = 0.8, Qb = 0.04 and finally Hq = 10^ h m/s/Mpc, with 
h = 0.72. The dark energy equation of state is set to w = -0.95 
and the sound speed is equal to the speed of light, Cs = c, such that 
there is no dark energy clustering. 



2 COSMOLOGY 

2.1 Dark energy cosmologies 

The dynamics of a spatially flat Friedmann-universe with dark mat- 
ter and dark energy is described by the Hubble function H(a) - 
d\na/dt, which is given by 

— 2- = — + (1 -QJexp 3 dlna 

with the matter density parameter and the dark energy equation 
of state function w(a). The value w = -\ corresponds to the cos- 



i{l +w{a))\. 



(1) 



mological constant A. Comoving distance x and scale factor a are 
related by 



da 



(2) 



which yields distances in unit of the Hubble distance;^// = c/Hq. 
For the galaxy redshift distribution n(z)dz, we use a standard shape 

n(.) = «o(|fexp(-(lf]d. w.th l = |r(|), (3) 

with p arameters zo = 0.64 and j3 = 3/2, as forecasted for EU- 
CLID jAmara & Refregier"2007S. The distribution can be rewritten 
in terms of comoving distance using the relation p(z)dz = pix)^X 
with dzldx = H{x)lc. 



2.2 CDM power spectrum 

The spectrum P{k) describes the fluctuation amplitude of the 
Gaussian, statistically homogeneous density field 6, {d{k)d{k')) - 
{2n)^6D{k + k')P{k), and is given by the ansatz 



(4) 



with the transfer function T{k). This transfer function is approxi- 
mated by 



T{q)- 



ln(l + 2.34^) 
2.34^ 



(1 + 3.89^ + (16.1^)2 + (5.46^)^ + (6.71^)^) 



(see'Bardeen et al."l986'). The wave vector k - q Y is measured in 
units of the shape parameter F. Sugiyama (Il995h describe correc- 
tions due to a non-zero baryon density Q^, 



F - Q,nh exp -Qi 



1 + 



^/2h 



(6) 



The spectrum P(k) is normalised to the variance erg on the scale 
R^S Mpc/h, 



2 _ r k^dk 



P(k)W\kR) 



(7) 



with a Fourier transformed spherical top hat filter function, W(x) = 
3 j](x)/ x. jfjx) is the spherical Besse l function of the first kind of 
order i jAbramowitz & Stegunll 19721) . 



2.3 Structure growth with clustering dark energy 

The growth of the density field in the linear regime, 6{x,a) - 



lows as a solution to the growth equation ( 


Turner &Whitelll997l: 


IWans & Steinhardtlll998l:lLinder & JenkinsI 


2003b. 



d^ ^ + i (3 + ^^^^ \ A/) (a) - — 
da^ ^ a\ d\na ) da ^ 2a^ 



2.4 Lagrangian perturbation theory 



Q^(a)D+(a). 



(8) 



The peculiar motion of galaxies can be described using Lagrangian 
perturbation theory (LPT) if the flow of dark matter and of the 
advected galaxies is irrotational and nonlinearities are weak. In 
this limit, galaxies follow straight lines given by the gradient 
of th e Zel'dovich potentia l Oi to first or d er (ILPT, Zel'Dovichl 
19701: lDoroshkevichlll970l : iBucherJ [19891 : iMoutarde et al...l99L : 
Bernardeau et al.ll2002h . 
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x^x-Z)i(a)VOi 



(9) 



where Oi is the solution to the Poisson equation, AOi = AO = 
6. This solution can be improved by adding second order cor- 
rections to Lagrangian perturbation th eory f 2LPT. iBuchertI 1 1 9941 : 
lMelottetalJll995l : iBouchet etal1l995h . 

X ^ X- Di(a)VOi + D2(a)V02 

with the second order potential 
iBouchet etaLlll995h . 



(10) 

$2 teuchert et alJ Il994 

(11) 



The time dependences are given by Di(a) = D+(a) and £>2(a) = 
-3/7Z)i(fl)fl~'^ '"'^ (for a low Pm- cosmology with a cosmological 
constant A, see IBouchet et The solution to the latter re- 

lation can be written down in Fourier- space, where the products of 
tidal fields become convolutions, 

^2 = -^ / 0pJ]Qij(k,k')6(k')6(k-n (12) 



where the mode coupling function Qij(k, k') becomes: 
^ _ {k'f^jk' - k)] - kjkjjk - k'Uk - k')j 

^ ^ " (k'm-ky • 

Spectra of the potentials and (^2 can be defined by 



with Pfik) = 



{2nfSo{k + k')P^^{k), 



1,2, 



(13) 



(14) 



P{k)/k'^ as a consequence of the Poisson equation 



and with P)p (k) which can be derived to follow 

\2 



^•'^'^^^l/^ l]Qdk\k-k') P(\k'\)P(\k-k'\)(l5) 



by application of the Wick- theorem (for a proof, see lDurreiil2008 
The integration is most efficiently carried out using cylindrical co- 
ordinates aligned with k such that d^k' = 2;r(^')^d^'dcos ^ using 
azimuthal symmetry, with 6 being the angle between k and k\ 

Fig- Ogives an impression of the spectrum P^^\k) and of the 
2LPT-corrections Pf(k) relative to ILPT. We plot k^P^^ik), / = 1, 2 
which is equal to the CDM spectrum P{k) for the ILPT result due 
to the Poisson equation. The 2LPT-spectrum is smaller on almost 
all scales by up to an order of magnitude and is only similar in 
amplitude on spatial scales of about 1 Mpc/h. 



2.5 Angular momentum from tidal shearing 

Angular momenta of dark matter haloes are introduced by tidal 
shearing, where the di ff'erential mo t ion of a protohalo gives rise 
to a torquing moment (Hovle"l949l: ISciamal[l955l : |Peebleslll969l : 
IBoroshkevich 1970; White 1984): 



3 dD+ 

La - a H(a)— €afiy / ^ I I 



(16) 



This relation reflects the interesting misalignment property of the 
shear and i nertia eigensystenis nece ssary for angular momentum 
generation (ISchafer&Merkelll201lh : Only the antisymmetric ten- 
sor = Tjdi^fis^sy ~ ^pshy)!^ is relevant for the angular momen- 
tum. La oc X^^, because the contraction of the symmetric tensor 
XJ^ = ^s(^j3s^sy + ^/sshy)/^ with the antisymmetric ea/ss vanishes. 
The antisymmetric tensor X_ is equal to the commutator [I^g, O^^] 
which suggests that for angular momentum generation, the tidal 




wave vector k [Mpc/h] 



Figure 1. CDM spectra k^P^^(k), i = 1,2 which are employed in La- 
grangian perturbation theory for displacing the galaxies, for ILPT (solid 
line) and 2LPT (dashed line). The ILPT result corresponds in this represen- 
tation to the CDM spectrum Pik): k'^P%\k) = P{k) 



shear and the inertia are not allowed to be simultaneously diago- 
nalisable and may not have a common eigensystem. 

Angular momenta L are described as being coupled to the tidal 
shear by means of a Gaussian random process p{LYbij)&L involv- 
ing tidal fiel ds O/, shaping the covariance cov(L)/y of the Gaussian 
distribution (iLee & Penll200lh . 



COV(L),-y = {LiLj) : 



I + a 



(17) 



with the misalignment parameter a, which describes the average 
orientation of the protohalo 's inertia to the tidal shear eigensystem. 
a has been measured in numerical simulation to be close to 0.25 
which we will assume in this work. O is the unit normalised trace- 
less tidal shear with the properties tr(0) = and tr(O^) = 1. This 
description is valid on scales where the correlations between tidal 
shears are negligible. 



2.6 Intrinsic ellipticity correlations 

Ellipticity correlations between galaxies are traced back to cor- 
related angular momenta of their host haloes. CDM haloes ac- 
quire their angular momentum by tidal shearing and due to the 
fact that neighbouring galaxies experience correlated tidal fields, 
their angular momenta are correlated in consequence. The angu- 
lar momentum L in turn determines the angle of inclination at 
which the galactic disk is viewed, and ultin i ately the ellipticity e 
(iHeayens et aDl200Ql: Erittenden et all 1200 iL I2OO2I : iMackev et aP 
l2002l : lHevmans & Heavensll2003h : 



with 



1 + 1^? 



= 2a 



LxLy 

1 +L?' 



(18) 



with the angular momentum direction L - LjL and the coordinate 
system being aligned with its z-axis being parallel to the line of 
sight. A rotation of the coordinate frame by (p causes the complex 
ellipticity to rotate twice as fast, e exp(2i(^)6, in accordance 
with the spin-2 property of the ellipticity field, or is a free parameter 
weakening the dependence between inclination angle and ellipticity 
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for thick galactic disks and has been determined to be or = 0.75 in 
the APM sample (Crittenden et al. 2001). 

In this work we use the angu lar momentum-based e lliptic- 
ity correlation model proposed by ICrittenden et al. 1 (l200lh . who 
trace ellipticity correlations back to tidal shear correlations using 
the conditional probability distribution p(L\(^ij)dL introduced by 
iLee & Pe n (2001): In this model, the distribution p(L\(^ij)dL is as- 
sumed as being Gaussian which is then being marginalised over the 
magnitude of the angular momentum vector, retaining only its di- 
rectional dependence. Writing down the ellipticity components as 
a function of the angular momentum direction and employing the 
covariance (LiLj) as a function of the squared tidal shear tensor, 
as advocated by Lee and Pen, it is possible to relate the tidal shear 
correlations to the spectrum of the density field. With this relation, 
one can write down a correlatio n function of the ellipti city field as 
a function of moments ^n(r) (see lCrittenden et aDl200lh of the tidal 
shear field and finally to carry out a Limber projection for obtaining 
the angular correlation function. For the parameter a we chose the 
value 0.25 supported by numerical simulations. 

Ellipticity correlations between two points 6i and 02 separated 
by the distance 6 are described in terms of two correlation functions 

UO) = {e{9,)e{92)) = (6,(^1)6,(^2)) + (6x(^i)6x(^2)> (19) 

= (6(^i)6(^2)> = (6,(^l)6,(^2)>-(6x(^l)6x(^2)> (20) 

which are formed from the variances of the ellipticity components 
6+ and 6x using {e+Sy) - 0. They can be transformed to the spectra 
C^(^) and C^(^) of the gradient and vorticity modes of the elliptic- 
ity field. 



cm 



n j ede[uo)Jom + UO)J4m], 



(21) 



(22) 



by Fourier transform (Kaiser 1992; Schneider et al.' l2002l : 
ISchneider &" Kilbinger 2007; Fu & Kilbinger 2010). Fig. [6] shows 
intrinsic ellipticity spectra C^(^) and C^(^) for the EUCLID galaxy 
sample with its median redshift at Zmed = 0.9. The spectra are con- 
stant and equal in amplitude up to multipoles of ^ ^ 100, indicating 
the absence of correlations such that on each scale on measures the 
variance of the uncorrected ellipticity field. Correlations become 
important on angular scales £ ^ 300 where the spectra level off" and 
decrease from multipoles of ^ ^ 3000 on very rapidly. In the peak 
region, the ellipticity £-modes have an amplitude larger than the 
5-modes by about an order of magnitude. 



3 EVOLUTION OF ELLIPTICITY CORRELATIONS 

3.1 Analogy between ellipticities and polarisation 

The evolution of the angular ellipticity spectra due to peculiar mo- 
tion of the galaxies are described in our model by drawing an anal- 
ogy to lensing of the polarisation modes of the cosmic microwave 
background. Both the galaxy ellipticities and the Stokes-parameters 
of the CMB -polarisation form a tensorial spin-2 field, which means 
that rotations of the coordinate frame by an angle (p give rise to a 
transformation of the tensor components as 6 ^ exp(2i^)6 and 
P exp(2i(^)P, when the ellipticity is written as a complex ellip- 
ticity 6 = 6++ i6x and the polarisation tensor P is composed of the 
Stokes parameters Q and U according to P = U + iQ. 

Peculiar motion as well as gravitational lensing introduces a 



shift in the position by an angle a such that the ellipticity e is not 
observed at the position 6 where the galaxy was formed, but rather 
has been displaced 6(6) e{6 + a). The correlation properties of 
such a distorted field can be computed using the formalism devel- 
oped for CMB lensing, which allows the computation of correlation 
of the lensed polarisation field, P{6) P(6 + a), where a refers 
now to the lensing deflection angle. Our formalism will be built in 
complete analogy and computes the shifting angle from the pecu- 
liar velocity, which in turn is derived from a velocity potential using 
Lagrangian perturbation theory for the description of peculiar mo- 
tion. 



3.2 Formalism for displacing the ellipticities 

By drawing analogies between the peculiar motion of galaxies 
causing displacements in the ellipticities, e(6) e(d + a) and 
the lensing of the polarisation of the CMB, P(6) P{6 + a) it 
becomes possible to derive spectra of the evolved ellipticity field. 
Peculiar motion by Z)+(a)VO changes the position of a galaxy by a 
shifting angle a - D+VO/;^ if the galaxy is situated at a comoving 
distance x- The angular displacement field a can be derived from 
a displacement potential ij/ - D+O/x^ by angular derivation, such 
that a = Voif/, because = x^- Generalising this argument to a 
galaxy population which is described by a normalised distribution 
n(x)dx in comoving distance x one obtains an expression for the 
angular displacement potential, 

n(x)D^ 



f 



dxW^(x)<b with W^(x)- 



(23) 



which replaces the lensing potential in the case of gravitational 
lensing of the CMB. The statistical properties of j/^, which is a 
Gaussian random field, are described by the spectrum C^{1), 



CM) : 



(x)P^(k = tlx) 



(24) 



which results from carrying out a Limber-projection of ij/. The spec- 
trum Ca{i) is related to C^{i) by Ca{i) - PCfi,(£) as a consequence 
of the relation a = Voif/. 

The angular spectrum C^^(£) of the displacement potential if/ 
resulting from the Limber-projection of Pq,(k) is depicted in Fig.|2] 
along with the spectrum Ca(0 = i'^C^p{i) of the displacement angle 
a. Clearly, the ILPT-result dominates over the 2LPT result by more 
than one order of magnitude, as already suggested by Fig. [T] The 
similarity of the plot to the analogous quanities in CMB -lensing is 
striking. 

Correlations between the components of the shiftin g angle a 
at two positions 6i and 62 are described by (ISeliaklll996h 

{ai{9,)aj{92)) = ^-CoiO) - €2(6) 6^6 (25) 

with 6 - 62 - Oi, and correlation functions of the displacement 
angle which are defined as 



rfid£ 

Co{e) = J -^c^(^)Jom 

and 

—cmJim- 

We introduce the abbreviation 0-^(6) = Co(0) - Co(6) in com- 
plete analogy to CMB -lensing for describing uncorrected displace- 
ments. The characteristic function of a Gaussian displacement field 
a would then be: 



(26) 



(27) 
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C^(l), 1LPT 

C (I), 1LPT 

_ _ C^(l), 2LPT 
C (I), 2LPT 




30 1 00 300 

multipole order £ 



Figure 2. Angular spectrum C^(^) (green line) of the displacement potential 
4f and the spectrum Ca(i) = fiCtf^if) (blue line) of the displacement field 
a = Veip, for ILPT (soHd line) and 2LPT (dashed line). 




10 10 10 

separation angle 6 [arcmin] 



Figure 3. Correlation functions (j'^iO) = Co(0) - Co(e) (blue line) and €2(6) 
(green line) as a function of separation angle 6, for ILPT (solid line) and 
2LPT (dashed line). 



<exp ae [a{6,) - a{62)])) = exp |- [-o-\e) + cos 2^,C2(^)jJ .(28) 

In the case of CMB-lensing, non-Gaussian contributions have been 
shown to have neg ligible e ffect on the deflection angle statistic 
(ICarbone etalj[2009l : IMerkel & Schafer 2011) and in the case of 
weak cosmic shear, analogous arguments about th e sparcity of 
strong deflections apply equally (iHamana et alJl2005h . 

Fig. [3] shows the quantities 0-^(6) = Co(0) - €0(6) and €2(6) 
used in this formalism, for both ILPT and 2LPT. Again, we would 
like to draw the reader's attention to the similarity between our re- 
sults and the formally equivalent result in CMB-lensing and to the 
domination of the ILPT results over the 2LPT spectra. 

The correlation properties of the shifted ellipticity field can be 
described using the two correlation functions ^+(6), 



= {€*(x + a)€(x' +a')) 
C(0) = (exp(-4i0^)6(x + a)6(x' + a')) 



(29) 
(30) 



where the points at which the ellipticities are observed, are shifted 
by exactly the angle a. Substituting the correlation function for the 
deflection angle in the Fourier- transforms of the above expressions 
yields the correlation functions ^+(0) of the shifted ellipticity field. 
They can be transformed to £^-mode and 5-mode spectra with the 
standard transformation written down in eqns. HTl and (|22] ). 

In summary, the £-mode and 5-mode spectra of the shifted 
ellipticity field can be written concisely in a matrix notation: 



W.(£, £') 



c|(r) 



(31) 



This notation shows explicitly the mixing between scales due to the 
convolution integral and the conversion between C^(0 and C^(0 
under the influence of W-{£, i'), which is the non-diagonal entry of 
the mixing matrix. The kernels W+{i, I' are given by 



J') = ^ J" Ode [Jo{£e)A{£\ 6) - J4(W)B(f, 6)] , 



with the functions 
A(i, 0) = exp 



B{i,e) = exp - 



2 



Jo(£,e) + -C2(e)um 



J4(ie) + -C2(e)Jsm 



(32) 
(33) 

(34) 
(35) 



which describe uncorrected shifting due to cr^(6) and correlated 
displacements due to €2(6). We abbreviated Js(x) = J2(x) + Jei^)- 
In the limit of no shifting, €0(6) = €2(6) = such that W+(^,f) = 
6(£ - and W-{i,i') = 0, due to the orthogonality relations of 
the cylindrical Bessel functions. 



Ode Jn{ie)Jn{i'e) = ^d^ii - O- 



(36) 



In this limit, the convolution is reduced to a Dirac (5z)-function and 
the mixing matrix is the unit matrix, so that the £^-mode and B- 
mode amplitudes are conserved. We have verified that higher-oder 
corrections arising in the transformation of correlation functions 
do have a negligible eflfect for the evolved e llipticity correlations 
dChallinor & Lewisll2005l : lLewis & Challinoi][2006) . 



3.3 E 1 5-mode conversion 

Figs, m and [5] show the mode coupling kernels W+{i,i') and 
W where for simplicity we focus on ILPT because the con- 
tributions due to 2LPT are comparatively small. From Fig.lUwe see 
that the power of the W+(^,^') -kernel is mainly distributed along 
the diagonal and increasing with multipole number, with maximum 
contribution from 300 ^ 3000. The oflf-diagonal contribution 
creates a convolution (eqn.[3TJ between the spectra at diflferent mul- 
tipoles, mediated by W+{i,i'). In contrast, the mode coupling ker- 
nel W-{i, i') (Fig. [5]), which is responsible for the E'/^-conversion, 
shows a lateral pattern which is three orders of magnitude smaller 
in amplitude and decreasing with higher multipole numbers {i, i'). 



3.4 Ellipticity spectra 

The final result is given in Fig. (6] which compares the initial ellip- 
ticity spectra C^(0 and C^(0 of the ellipticity field as predicted by 
the angular momentum model, and the evolved spectra C'^{€) and 
C'^{€) due to peculiar motion. For comparison with weak lensing. 
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Figure 4. Mode coupling kernel (ii') x W+(i,i') used in the transforma- 
tion of the ellipticity spectra. For simpHcity, we focus on ILPT because the 
contributions due to 2LPT are small. 




100 300 
multipole I 



1000 3000 



Figure 5. Mode coupling kernel i£f)xW-i£, i') responsible for the E ^ B- 
mode conversion in the ellipticity field. Again, we show the results for ILPT 
because the contributions from 2LPT are small. 



we plot the weak convergence spectrum Q(^) expected from the 
EUCLID galaxy sample in compariso n, for a nonlinear C DM spec- 
trum (using the parameterisation bv ISmith et alJ 120031) . The first 
observation is that ellipticity correlations reach amplitudes similat 
to those of the weak lensing convergence in the nonlinear part cor- 
responding to amplitudes ^ ^ 300, and that the intrinsic £-mode 
spectrum C^(0 is larger than the 5-mode spectrum Cy{i) by about 
an order of magnitude in this regime. On larger angular scales, there 
are no appreciable ellipticity correlations and one efifectively ob- 
serves the variance of the ellipticity field for uncorrected objects. 
Consequently, the spectra have identical amplitudes and are eflfec- 
tively constant. In this regime, the shifting eflfect is not able to af- 
fect the galaxies, which is a well-kno wn result in CMB -lensing , 
where scale free- spectra are invariant (iLewis & Challinojl200d) : 
The mode-conversion mechanism is uneflTective if the spectra are 
equal, C^(0 = and the convolution with W+{i, i') is not able 

to redistribute amplitudes. In contrast, both spectra are aflfected on 
multipoles i > 1000, where in particular C'^{i) has decreased rela- 
tive to C%{i). 

Fig. |7] compares the relative magnitude of all spectra as a 
function of multipole i. The plot shows the relative ratio of the 
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Figure 6. Ellipticity spectra C|(^) (blue line) and C^(^) (green line) as 
predicted by the angular momentum model with a = 0.25 and the disk 
thickness parameter set to a = I (dashed line), and the evolved ellipticity 
spectra (solid line) where the displacements were computed by ILPT. For 
comparison, we plot the spectrum Q(^) of the weak lensing convergence for 
a linear (black dashed line) and nonlinear (black solid line) CDM spectrum. 
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Figure 7. Ratios C^(^)/C^(^) (blue sohd line), C^(^)/C^(^) (green sohd 
line), Cj^{i)IC^ii) (blue dashed line) and C]^{i)IClii) (green dashed line) 
with all displacements following from ILPT. 



evolved and initial £^-mode and 5-mode spectra. As already indi- 
cated by Fig. [6] we see a significant decrease for / > 1000 of up 
to 10% for the E- and 60% for the 5-modes at / ^ 3000. The ra- 
tios C^(0/C^(0 and Cl{£)IC'^{£) of intrinsic and evolved spectra 
are similar up to multipoles of ^ ^ 1000, where they separate and 
indicate that the newly generated 5-modes are small and that the 
5-mode spectra are more strongly aflfected. For EUCLID 's weak 
lensing application, changes in the ellipticity spectra are aflfecting 
scales where the shape noise starts dominating, but for shallower 
surveys, lower multipoles would be aflfected by the peculiar motion 
eflfect. 

Finally, Fig. [8] gives an impression of the mode conversion 
mechanism, where we plot evolved spectra C'^{i) and C'^{i), when 
the £-mode or the 5-mode in the initial spectra was deliberately 



© 2008 RAS, MNRAS 000,[T}i] 



Evolution of ellipticity correlations 7 



10"^ 




multipoie order i 

Figure 8. Contributions to the evolved ellipticity spectra C'^ii) (blue lines) 
and C'^{i) (green lines): no initial 5-mode spectrum, C^(^) = (dashed 
lines) and no initial £^-mode spectrum, C^(^) = (solid lines). 

set to zero, i.e. C^(0 = in the first and C^(0 = in the second 
case. Even in the absence of a particular initial mode we observe 
power in the corresponding evolved spectrum, as a consequence of 
£'/5-coupling introduced by peculiar motion. 



for C^(^) and 60% for C^(^). The mode-conversion mechanism is 
comparatively weak and we tested it by deliberately setting the ini- 
tial spectra C^(0 and C%{i) to zero. 

(iii) Second order corrections in the dynamical model were 
found to be negligibly small in comparison to first order Lagrangian 
perturbation theory. Likewise, we made sure that higher-order cor- 
rections in the transformation of the ellipticity spectra had a minor 
eff'ect on the evolved ellipticity spectra. 

We conclude that in principle the dispersing eff'ect of peculiar mo- 
tion weakens intrinsic ellipticity correlations and make them less 
troublesome for analysing weak lensing data. For the case of EU- 
CLID we see changes in the spectra on scales where the shape 
noise is already dominating. A natural extention to this investi- 
gation would comprise the shifting and distorting eff'ect of weak 
gravitational lensing, and ultimately the usage of analysis methods 
conceived for the polarisation of the CMB for investigating intrin- 
sic ellipticity correlations. 
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4 SUMMARY 

The topic of this paper is the evolution of intrinsic ellipticity corre- 
lation between galaxies due to peculiar motion. Intrinsic ellipticity 
correlations are derived in the framework of angular momentum 
models, which explain these correlations by correlated tidal shears 
experienced by the protohaloes in acquiring their angular momenta. 
Because the symmetry axis of the galactic disk is related to the an- 
gular momentum direction of the host halo, correlated angular mo- 
menta give rise to correlated angles of inclination and hence corre- 
lated ellipticities. 

(i) Peculiar motion of galaxies changes the correlation proper- 
ties of the ellipticity field by displacing the galaxies and distorting 
the ellipticity field. We describe the peculiar motion by Lagrangian 
perturbation theory and derive corresponding displacement angles 
along with their statistical properties for the EUCLID galaxy sam- 
ple. The formalism for evolving the ellipticity spectra uses an anal- 
ogy to the formalism describing lensing of the CMB polarisation 
spectra. Both quantities, the ellipticity field as well as the polarisa- 
tion field, have the same symmetry properties, being of spin 2. The 
loci at which ellipticities and polarisations are measured are dis- 
placed by peculiar motion in the first and by gravitational lensing 
in the second case. Because the peculiar motion field in the quasi- 
linear regime is a flow resulting from a velocity potential which 
corresponds to the lensing potential, is it possible to derive all nec- 
essary quanities in complete analogy. 

(ii) Peculiar motion has two eff'ects on the ellipticity spectra: 
There is a convolution of the spectra and a conversion between E- 
modes and 5-modes of the ellipticity field. Both eff'ects become im- 
portant on angular scales £ > 1000, because on smaller multipoles, 
the spectra are eff'ectively constant and equally large. In particular 
the spectrum C^(^) is strongly aff'ected and looses amplitude: For 
the EUCLID galaxy sample we measure decrements by about 10% 
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